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Abstract 

We estimate the amplitude of thermal fluctuations by calculating the typical size of subcritical bubbles in cosmo- 
logical electroweak phase transition and show that this thermal fluctuation effect drastically changes dynamics of the 
phase transition from the ordinary flrst order type with supercooling. From this fact, we conclude that the standard 
electroweak baryogenesis scenario associated with such a flrst order transition does not work in the minimal standard 
model in certain conditions. 



Recently electroweak baryogenesis scenario has been actively investigated [1] [2]. In this scenario the phase transition 
must be of flrst order with supercooling in order to attain a non-equilibrium state for baryon number non-conserving 
processes. In one-loop appoxmation, the flnite temperature effective potential always yields a barrier between the 
false and the true vacua, i.e., a flrst order phase transition [3]. However, if thermal fluctuation is too large, the 
perturbation scheme breaks down and its prediction of a flrst-order phase transition becomes suspect. Although the 
dynamics of the phase transition is of primary importance in such a cosmological context, the analyses so far done 
have mainly been conflned in other aspects, namely, spontaneous baryon-number generation in thick walls [1] and 
charge transport mechanism in thin walls [2]. In fact, if the transition is not of flrst order with supercooling, then 
these mechanisms do not work, invalidating the entire scenario. In the present paper we consider effects of thermal 
fluctuations on the electroweak phase transition paying attention to the role of subcritical bubbles. 

So far Gleiser et al. [4] and Dine et al. [5] have estimated the amplitude of thermal fluctuations and obtained 
the result that they are dominant if > 80GeV [4] and therefore the critical bubble cannot be borne due to 
inhomogeneities of background flelds. As the experimental constraint is niH > 60GeV [6], the parameter range 
allowing electroweak baryogenesis is still open. However, in their analyses, they assumed that the spatial scale of the 
thermal fluctuation is equal to the correlation length, that is, the inverse mass scale. In the microscopic level, the true 
nature of the dominant thermal fluctuation would be a perpetual creation and annihilation of spherical subcritical 
bubbles. Thus one should identify the typical size of the bubbles with that estimated by the statistical ensemble 
averaging. What we would like to claim here is that the amplitude of the thermal fluctuation sensitively depends on 
the size of the subcritical bubble. A small change of the size results in a drastic change of the phase transition. 

To derive the characteristic scale of a subcritical bubble we shall construct effective Lagrangian for the size of the 
bubble R(t). In high temperature phase, the motion for the wall of asymmetric subcritical bubble is bounded from the 
above and then one can obtain the flnite thermal average scale R of R(t). On the other hand, in the low temperature 
phase, symmetric subcritical bubble is bounded from the above also and the size is almost the same as the above R. 
The correlation length is the maximum scale of thermal correlation and dose not necessarily coincide with the size of 
the subcritical bubble in general. This means that R is smaller than the correlation length and that the amplitude of 
thermal fluctuation will increase, invalidating perturbation analysis. Intuitively, this means that one cannot flnd the 
distinction between symmetric vacuum and asymmetric vacuum and that the spontaneous critical bubble nucleation 
out of the homogeneous symmetric state, which is usually assumed, cannot happen. 

In this paper we concentrate on the minimal standard model. The effective potential with one-loop and flnite 
temperature correction [3] is given by 
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and 

where we have taken niw = 80.6GeV, mz = 91.2GeV, mj = 174GeV and a = 246GeV [7]. In this case T2 = 
\/im'jj — 8Ba'-^)/4:D ~ 92.69GeV and the critical temperature at which the hight of the potential at two local minima. 



and <j>+ := ^[l + Jl - §^S^(T^ - Tf)], degenerates is = Ts/^l - E^/\tD ~ 93.43GeV, where 



2At ' V 9E^T 



5 = ^ |6(-^r + 3(^r - 12(^)^1 ~ -0.00456. (5) 



The temperature at which the global minimum appears is Ti = ^2/^1 - 9E'^/8XtD ~ 93.52GeV. 

We have also chosen the minimal value niH = 60GeV for Higgs mass permitted from experimental constraint. As 
the thermal fluctuation is an increasing function of the Higgs mass, it is sufficient for our purpose if thermal fluctuation 
is shown to be too large with niH = 60GeV. On the other hand, in order to avoid baryon number destruction due to 
sphaleron-mediated processes after the phase transition, we need < 55GeV in one-loop perturbation theory [8]. 
Further, Shaposhnikov has obtained the result niH < lOOGeV in lattice calculation [9]. Recently, the loop correction 
for the sphaleron transition rate have been estimated and in this case niH < 63GeV [10]. Thus, the possibility for the 
explanation of the present net baryon-number remains at present time. 

As we have already mentioned, the typical size of thermal fluctuation is a calculable quantity. Since one cannot use 
the instanton method for them, we must invent new method. The subcritical bubble which minimizes the free energy 
would have the 0(3) symmetry. Hence in the Lagrangian, 

L= [d'jU^-hvct^r-VTW], (6) 
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we adopt the following ansatz for the bubble conflguration 



(f) = (f)+exp ^--j^^^, r:=|x|, (7) 

where R(t) is the size of subcritical bubble and is time dependent. Without any deflnite argument, previous authors 
have chosen the correlation length i{T) = 1/^2D{T'-' - T|) as the typical size R{t). However, in the microscopic 
level, the bubble has dynamics. Inserting (7) into the above Lagrangian (6) yields an effective Lagrangian for R(t): 

157r3/2 . '^7r3/2 tt^/^ 7r3/2 tt^/^ 

" 16^2 + 4^2 + 2V2 + 3^3 + 32 

Applying the variational principle on R(t) we obtain the equation of motion for R(t): 

^ + 7^(^)2 + ?4 + t^iT^ - Ti)R - ^ETct,+ R+^ct,lXTR = 0. (9) 
dt^ 2R^ dt ' b R b ^ ^' 15^3 10^2 "^ 

In the vicinity of the turning point \dRl dt\ <^ 1, one can omit the second term of the left-hand-side and then we get 
an energy equation by multiplying dR/dt and integrating over t, which reads 
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where Eg is a constant of motion and Rc := ^ Eg /U , T) is the turning radius. The potential U{R,4',T) is 
a monotonicaUy increasing function of R and therefore the wall motion is restricted from the above in the high 
temperature phase(T > T^). The existence of the upper limit will give the typical size of the subcritical bubble. 

What we would like to obtain is the typical size of a subcritical bubble. In order to estimate it we must construct 
the Hamiltonian which describes the dynamics for R{t). The canonical momentum is defined by 
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The effective Hamiltonian then becomes 
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Thus the thermal average of R is given by 
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(14) 



As a ~ 0.015 ^ 1 and a / [3 ~ 0.5, the last appoxmation is certainly justified. Further, one must treat carefully the 
measure. As 
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and 
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the original measure becomes 



dRdP / rf-^a; 



d((t),P4,) 

''d{R,p) 



dRdP. 



(17) 



Thus, naively used measure in eq.(14) is justified. Moreover, since the three-dimensional volume integral is done in eq. 
(17), weight associated with three volume has already been taken into account. The temperature dependence of eq.(14) 
is depicted in Fig.l. As one can easily see, the above result gives an apper bound for R. In fact, at T = and T = Ti, 
we numerically obtain the more accurate result that R{T = T^) = 0.012GeV~^ and R{T = Ti) = 0.0202GeV~^, 
respectively. In Fig. 1, we can easily find R{T) < £(T) as expected. 

Next, let us estimate the amplitude of the thermal fiuctuation on the scale R, adopting the trial function 



3 



i^exp[-=2]- (18) 
R 



Then the free energy becomes 



rS7r3/2_ 7r3/2 on 7r3/2 _^ _3/2 

i^o(./'^,T) = [:L_R+_D(T^-Ti)R']ct,l - —ET^\r' + —\t^\r\ (19) 



Thus the RMS amplitude of 4> at the symmetric vacuum is 



Jd(f)A(PAe T 



(20) 



Its temperature dependence is depicted in Fig. 2. In the same way as in the case of R, the above approximation 
result gives the lower bound for y^^- At T = T,, and T = Ti, one obtains numerically ^J~<p(T = T^) = 36.2GeV 

and \J~<p{T = Ti) = 28.2GeV, respectively. This exceeds the first refiection point <j)^, = — \J ^^T^ — ^^''sAj "^^ ; 

which implies that the perturbation theory breaks down and that one can no longer conclude that electroweak phase 
transition is of first order. For the comparison we calculate the RMS by the same treatment as Dine et al [5]. According 
to their analysis, it becomes 

/^2\ T <i>l_ ^21) 
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and then \/J(f)^ ~ 4'+l'i- This value is slightly smaller than eq. (20), but exceeds in the most range of temperature 
as one can see in Fig. 2. 

To obtain the general view of the phase transition we must examine the dynamics of the symmetric sub critical 
bubbles. In the low temperature, symmetric bubbles are also always bounded, and this gives the typical scale of the 
subcritical bubble. We adopt the trial function for a symmetric bubble in the asymmetric phase as follows: 



^+[1--P[-^]]. (22) 
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The same procedure as the above discussion yields the energy equation near the turning point R, 

l(dR\'^ 
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and is the constant of motion. Note that the above energy equation holds in the case in which turning point 
exists. In the low temperature phase(T < T^), the potential U{R,4>+,T) is a monotonically increasing function of R. 
To estimate the typical size we construct the effective Hamiltonian for the variable R{t) and its canonically conjugate 
momentum P{t): 

We obtain the same result 
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^ ^^m^ (26) 

as the asymmetric bubble in the same approximation level. By the same argument with the asymmetric bubbles, this 
gives an upper bound for numerical results. 

Let us clarify the total dynamics of electroweak phase transition. The following argument is based on the work in 
ref. [11]. Although it is a rough estimate, the results give us qualitative picture which is necessary for understanding 
the phase transition. Then transition rate becomes 

r[o.,] = m(T)exp[-^], (27) 



and 



r[+^o]=m(T)exp[-^^], (28) 

with m{T) = ^2D{T'^ - 7|), where 

i^o(^+) = \t+^, \mT' - Tim' - + 71 ^^^^1 (29) 

and 

-7r3/2(l--L- JL+ 1 )Ay,/,-2]T3. (30) 
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For the more detailed analysis of phase transition the above formulae are not sufficient. In particular, note that the 
prefactor m(T) is determined from the dimensional argument although this term is important when we determine the 
freeze out temperature Tfo- The master equations for the number density of the symmetric (asymmetric) state, Nq 
(#+), are 



dNo 
dt 



~dt 



-ro#o + r+#+, (31) 



-r+#+ + ro#o, (32) 



respectively. The behaviour of the both transition rates can be seen in Fig. 3, where we can find that Tjo is barely 
larger than T2. However, this results is mainly due to the prefactor m(T). As the both transition rates are much larger 
than H = d/a during the phase transition, Tfo is probably lower than T2 if one can calculate the exact transition 
rate. For example, if we substitute m(T) by T, Tfg becomes lower than T'j. Thus the fraction of false vacuum for 
true vacuum Y = Nq / has the equilibrium value at T = T2 : 

^-^±-exp L^^o 



y = y''' = ^ = ^ = -^v\-^T^\^ (33) 



Its temperature dependence is seen in Fig. 4. When the temperature becomes Ti, the asymmetric vacuum appears. 
The fraction which the asymmetric vacuum occupies will be almost zero at this moment. The transition rate from 
asymmetric to symmetric vacuum exceeds during a small temperature difference from T = Ti. When the temperature 
becomes the critical temperature, the fraction Y becomes unity. If the freeze-out temperature is higher than T2, 
after the temperature becomes the freeze-out temperature the symmetric region will collapse by surface tension. The 
phase transition cannot be accomplished by critical bubbles. If the freeze out temperature is lower than T2, the phase 
transition becomes of the second order type. Here we remark that the last attempt exists and some new aspect have 
been obtained in the case with light Higgs mass(miy < 53GeV) by considering the kinematics of the subcritical bubble 
[12]. For such a case, destruction of symmetric bubble by thermal noise becomes dominant. 
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In this paper, we investigated the strength of the thermal fluctuation and the structure of the phase transition. 
Our main focus is the estimation of the size of sub-critical bubbles. According to our calculation, the typical size is 
smaller than the correlation length and the amplitude of thermal fluctuation increases compared with the previous 
estimation. For experimentally allowed values of Higgs mass, the amplitude of thermal fluctuation always exceeds the 
first refiection point and then perturbation scheme is not guaranteed. At the critical temperature, the fraction which 
symmetric and asymmetric vacua occupy is almost the same and therefore the critical bubble cannot be nucleated. 
Note that Borrill and Gleiser have obtained the results that the two phase mixing is completely attained in numerical 
experiment [13] and our discussion might be supported by this. 

When we obtain effective Lagrangian for the size R(t) we used the thick wall ansatz. In this case, one can easily 
confirm that the typical size R of subcritical bubble is comparable with the thickness of the wall. Thus, the thin wall 
is not suitable for subcritical bubbles. 

In the microscopic level, the subcritical bubble expands and then collapses and any departure from thermal equilib- 
rium is not provided. Thus, the baryon number cannot be generated in minimal standard model. We conclude that 
the baryogenesis is quite difficult in the minimal standard model, even apart from the fact that CP-violation source 
is in Kobayashi-Maskawa phase only which is quite small. Thus, we need to investigate "beyond the standard model" 
and this work will be reported in the near future. 
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Figure Captions 

• Fig.l: The thermal correlation length £ and the typical size R of subcritical bubbles. The unit of the vertical 
axis is GeV""'^ . 

• Fig. 2: The RMS of thermal fiuctuations. The unit of the the vertical axis is GeV. 

• Fig. 3: Transition rates. Both rates degenerate in this figure. 

• Fig. 4: The fraction of the symmetric phase. 
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